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ABSTRACT 


To investigate the role of wind forcing, bottom topography and thermohaline 
gradients on classical as well as unique features in the northern Canary Current system 
(NCCS), four experiments are conducted with a sigma coordinate primitive equation 
model. The first experiment, which investigates the pressure gradient force error, shows 
that velocity errors inherent in three dimensional sigma coordinate models can be 
successfully reduced from ~1 m/s to less than 0.5 cm/s in the NCCS. The second 
experiment, which investigates the effect of annual wind forcing on a flat bottom, 
accurately portrays classical eastern boundary current features as well as unique NCCS 
features associated with a large embayment (i.e., the Gulf of Cadiz), poleward spreading 
of Mediterranean Outflow, and the generation of Meddies. The additional effect of 
bottom topography in Experiment 3 shows that topography plays important roles in 
intensifying and trapping the equatorward current near the coast, in weakening the 
subsurface poleward current and in intensifying eddies off the capes of Iberia. The use of 
full instead of horizontally averaged thermohaline gradients in Experiment 4 highlights 
the development of the Iberian Current off the Portugal west coast, a feature not seen in 
the previous experiments. This shows that thermohaline gradients play an important role 
for the formation of the Iberian Current. 
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I. INTRODUCTION 


The Canary Current System (CCS) on the eastern boundary of the central North 
Atlantic is a classical eastern boundary current (EBC) system. Stretching from ~10°N to 
~45°N along the coasts of northwest Africa and the Iberian Peninsula (IP), it marks the 
closing eastern boundary of the North Atlantic Gyre. Typical of other EBCs, the mean 
equatorward Canary Current (CC) is a broad (-1000 km), relatively slow (-10-30 cm/s), 
yearlong surface flow extending to depths of -500 m (Wooster et al., 1976). The portion 
of the CC that stretches along the coast of the IP is often referred to as the Portugal 
Current (e.g., Tomczak and Godfrey, 1994). 

As in other EBCs, a poleward undercurrent exists near the coast beneath the CC 
(e.g., Meincke et al., 1975; Fiuza, 1980) as a relatively narrow (-10-40 km) and weak 
(-2-10 cm/s) flow, which is strongest between -100 and -600 m depth. The depth and 
strength of the undercurrent varies seasonally and latitudinally. In winter it shoals to the 
north near Cabo da Roca (see Figure 1 for geographic locations, and Figure 2a for 
bathymetric contours and coastline geometry for the region), and forms a third flow 
component commonly referred to as the Iberian Current (IC) (Haynes and Barton, 1990). 
A nanrow (25-40 km), relatively weak (-20-30 cm/s), seasonal surface current, the IC is 
found trapped near the coast against the shelf break (Fiuza, 1980; Frouin et al., 1990; 
Haynes and Barton, 1990) and can occasionally be seen as far south as Cabo de Sao 
Vincente (Batteen et al., 2000). 

A unique feature that distinguishes the CCS from other EBCs is the existence of 
Mediterranean Outflow (MO) through the Strait of Gibraltar into the adjacent Gulf of 
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Cadiz. A large embayment, the Gulf of Cadiz’s pronounced east-west coastline 
orientation results in weaker upwelling in the Gulf of Cadiz than to the north or south of 
the Gulf of Cadiz, due to the dominant equatorward trade wind direction. The Gulf of 
Cadiz also creates a large separation between the two west coast upwelling regimes so 
that no continuous flow between the two appears to exist (Barton, 1998). 

The salty MO plume that exits the Mediterranean Sea through the Strait of 
Gibraltar is diluted, thickens, and becomes vertically differentiated into two distinct cores 
as it flows westward into the Gulf of Cadiz (Iorga and Lozier, 1999). At ~7°W in the 
Gulf of Cadiz, both a shallow core at depths of -600-900 m and a deeper core at -1100- 
1200 m exist. Both cores continue to flow westward along the southern coast of Spain 
and turn poleward around Cabo de Sao Vincente (Amber and Howe, 1979; Iorga and 
Lozier, 1999). A third, shallower, poleward core of Mediterranean water has also been 
traced from the Strait of Gibraltar northward to ~38.5°N off western Portugal (Amber, 
1982). In addition, a climatological cyclonic circulation in the southwestern Gulf of 
Cadiz acts to spread salty MO south of ~34°N (Iorga and Lozier, 1999). 

The northern CCS (NCCS) is influenced predominantly by equatorward, 

upwelling favorable winds produced by the eastern half of the Azores High. Located in 

the northeastern part of the Atlantic Ocean, the Azores High is a semi-permanent 

subtropical high pressure system similar in nature and behavior to the North Pacific 

Subtropical High (Nelson, 1977). As such, the center of the Azores High migrates 

meridionally with the seasons, reaching its southernmost extent near 27°N in March and 

ridging north to ~33°N by August. The building and migration of the Azores High cause 

wind stress values over the NCCS region to vary temporally, resulting in both alongshore 
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and cross-shore variability of the predominantly equatorward winds along the west coast 
of the IP and the northwest coast of Africa. The east-west pressure contrast between 
Portugal and the center of the Azores High during summer is ~8 mb. In winter, this 
pressure gradient weakens to ~1 mb. As a result of this seasonally changing pressure 
gradient, considerably stronger northerly and northwesterly winds occur in the NCCS 
region in summer, while northerly and northwesterly winds become weaker to even slight 
southerly off the Iberian Peninsula in winter (Batteen et al., 2000). The shift in 
maximum wind stress also causes upwelling favorable winds to shift from ~27°N near 
the Canary Islands in January, to ~43°N off Portugal by July (Fiuza, 1982). 

Like other classical EBCs, observations of the sea surface in the NCCS region 
have shown highly energetic mesoscale features such as jet-like surface currents, 
meanders, eddies and filaments over the broad climatological mean flow of the CC. 
Satellite sea surface images have shown nearshore upwelling during periods of upwelling 
favorable winds with several narrow filaments of cooler water extending off the coast of 
the Iberian Peninsula (Fiuza and Sousa, 1989) and Cape Ghir in northwest Africa (Van 
Camp et al., 1991; Hagen et al, 1996). In these images, upwelling filaments often 
extend -80-150 km offshore, with alongshore spacing of -80-100 km between filaments 
(Haynes et al., 1993), and often terminate with dipole eddy pairs (Fiuza et al., 1982; 
Barton, 1998). Observations have also shown anticyclonic and cyclonic pairs of 
mesoscale eddies on the order of 100 km off the IP coast (Fiuza, 1984; Stammer et al., 
1991). These mesoscale features have been observed during periods of predominantly 
upwelling favorable winds and appear to be located near prominent coastline 

irregularities such as capes. These observations provide evidence that wind forcing along 
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with coastline irregularities appear to be important mechanisms in the formation and 
sustainment of many of the mesoscale features found in the NCCS domain and other 
EBC regions (Batteen et al ., 2000). 

Unique to the NCCS is the generation of anticyclonic submesoscale coherent 
vortices (SCVs) or Meddles. Numerical studies suggest that baroclinic instability of the 
northward dense plume of salty MO along the IP continental slope leads to the generation 
of Meddies (Kase et al., 1989). As a result of numerous observations over the past 
decade, the primary generation region of Meddies is widely accepted to be near Cabo de 
Sao Vincente, off southwest Portugal. Several different trajectories of Meddies have 
been observed, including a southwestward movement into the Canary Basin, and 
westward translations south of the Azores (Richardson and Tychensky, 1998). 

Over the last two decades numerous modeling studies have focused on the driving 
mechanisms of complex mesoscale activities in EBC regions, including upwelling 
filaments, highly energetic eddies, and meandering jets. For the California Current 
System, the classical EBC system, baroclinic instability, wind forcing, and coastline 
irregularities have been investigated as possible generation mechanisms (e.g., Ikeda et al., 
1984a, b; Batteen, 1997). 

For the NCCS, Batteen et al., (2000) used the results of four numerical 
experiments of increasing complexity to investigate the classical and unique features of 
the NCCS. All four experiments were run on a 3-plane. Experiment 1 included seasonal 
wind forcing on a straight coast, while Experiments 2 and 3 added the effects of irregular 
coastline geometry to the seasonal wind forcing. The most complex experiment. 
Experiment 4, added the effects of both thermohaline gradients and annual MO to the 
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seasonal wind forcing and irregular coastline geometry. Consistent with previous studies, 
wind forcing was shown to be a key generative mechanism for EBC mesoscale features, 
while capes were shown to be key areas for enhanced upwelling, extensive filaments, 
maximum current velocities and enhanced growth of cyclonic meanders and eddies. 
Unique features of the NCCS, including the development of anticyclonic meanders and 
eddies in an embayment like the Gulf of Cadiz and the generation of Meddies from MO, 
were also highlighted. 

Another study by Batteen and Murray (2000) investigated the role of planetary 
beta on classical as well as unique features in the NCCS. Four numerical experiments 
were conducted with varying Coriolis parametrizations (f-plane or beta-plane). The first 
two experiments used a closed boundary and annual salinity forcing for the MO. The 
latter two experiments used an open Mediterranean Sea at the Strait of Gibraltar and 
seasonal forcing for MO to permit a more accurate investigation of the subsurface 
spreading of MO and Meddies. All four experiments used seasonal climatological winds 
and seasonal thermohaline gradients along the western boundary to force the model. 
Experiments run on a 3-plane (Experiments 2 and 4) accurately portrayed classical 
eastern boundary current (EBC) mesoscale features. In addition, these experiments 
depicted unique NCCS features associated with a large embayment (the Gulf of Cadiz), 
poleward spreading of MO, and the generation of Meddies. Experiments run on an /- 
plane (Experiments 1 and 3) showed the unrealistic dominance of a continuously 
strengthening equatorward jet that inhibited the development of classical EBC and unique 
NCCS features. The complex upper layer and subsurface flow regimes of Experiment 4 
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most realistically portrayed currents, mesoscale features and Meddies similar to NCCS 
observations. 

The objective of this study is to build on previous studies by incorporating the 
effects of bottom topography . The Princeton Ocean Model (POM), a bottom following 
sigma coordinate model, was chosen for this study because it has been widely used to 
simulate coastal processes associated with continental shelf flows and bottom boundary 
layer dynamics. The results of several numerical experiments (see Table 1) are explored. 
Each experiment includes the effects of thermohaline forcing at all three open 
boundaries, and uses a beta plane. In Experiment 1 velocity errors produced by the 
pressure gradient force error, an error inherent in all three-dimensional sigma coordinate 
models, are investigated using the horizontally averaged climatology with bottom 
topography and no wind forcing. In Experiment 2 the horizontally averaged annual 
climatology is used with annual wind forcing on a flat bottom. Experiment 3 is the same 
as Experiment 2 except that bottom topography has been incorporated. To explore the 
role of bottom topography on the NCCS, the results of Experiment 3 are compared with 
the results of Experiment 2. Experiment 4 is the same as Experiment 3 except that full 
annual climatology is used instead of the horizontally averaged annual climatology. To 
determine the role of the full climatology, results of Experiment 4 are compared with the 
results of Experiment 3. 

This study is organized as follows. In section 2 we describe the numerical model 
and the specific experimental conditions. The results of the numerical experiments are 
presented in section 3. A summary is presented in section 4. 
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II MODEL DESCRIPTION 


A. DATA SETS 

The topographic data were obtained from the Institute of Geophysics and 
Planetary Physics, University of California San Diego (Sandwell D.T. and W.F. Smith, 
1996). The data set has a resolution of 2 minutes (e.g., 1/30 of a degree) and is a 
compilation of 30 years of data soundings obtained by ships. Where the ships’ data is 
sparse, altimetry information was used to interpolate soundings . 

Annual temperature and salinity values were obtained from Levitus and Boyer 
(1994) and Levitus et al. (1994). The data uses a 1 by 1 degree horizontal resolution at 
the vertical levels shown in Table 2. 

For wind forcing, climatological wind fields were obtained from the European 
Centre for Medium Range Weather Forecasts (ECMWF) near-surface wind analyses 
(Trenberth et al., 1990). The data uses a 2.5 by 2.5 degree grid. 


B. PRE-PROCESSING 


The original topography (Figure 2a) was interpolated with a two-dimensional 
(2D) linear interpolation filter to the resolution used in the POM model, i. e., 3 by 3.7 km 
near the coast and 6 by 7.4 km away from the coast (Figure 2b), with a total of 287 by 
241 points. The highest resolution was used where the values of the ‘slope parameter’ 


(defined by Mellor, 1998, as 


1*1 


2 *H 


where H is the average depth and SH is the 


difference in depth between two adjacent cells), were the largest in both the latitude and 


7 



longitude directions. Since over much of the topography the slope parameter was larger 
than 0.2, which is the suggested maximum value to be used in sigma coordinate models 
(Mellor, 1998), the topography had to be smoothed with a linear 2D low-pass filter in 
order to meet this criterium. The new depth of each point calculated with this filter was a 
non-weighted average of 15 by 15 points surrounding the point. Subsequently depths 
greater than 2500 m were reassigned to depths of 2500 m, land was assigned the depth of 
10 m (to avoid divisions by zero in the model) and the Strait of Gibraltar was closed. The 
new topography is shown in Figure 2c. 

The annual temperature and salinity values were interpolated for the horizontal 
spatial resolution of the model and for the 21 vertical sigma levels (Table 3 and Figure 3) 
with a three-dimensional (3D) linear interpolation scheme. This had to be done separately 
for smoothed topography and for flat bottom due to the change in vertical levels for flat 
bottom and for topography. Temperature fields at sigma level one and several cross- 
sections are shown in Figures 4a to 4g. 

The daily seasonal winds were averaged over time in order to obtain the annual 
non-weighted average wind vector field (Figure 5). The wind vectors were interpolated 
for the horizontal spatial resolution of the model with a 2D linear interpolation scheme. 
The components of the wind stress were then calculated. 

C. BRIEF MODEL DESCRIPTION 

The Princeton Ocean Model, POM, a well documented model (e.g., Blumberg and 
Mellor, 1987; Mellor, 1996), was used in the model studies. POM is a primitive equation, 
free surface model with a second-moment turbulence closure scheme (Mellor and 
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Yamada, 1982) that, through the use of bottom-following sigma levels, can realistically 
simulate processes associated with continental shelf flows and bottom boundary layer 
dynamics in local domains (e.g, bays, estuaries and coastal regions). Recently, the model 
has been used sucessfully to simulate decadal processes in entire ocean basins (see Ezer 
and Mellor, 1994,1997). 

As described earlier, the resolution of the horizontal orthogonal grid varies 
between 3 by 3.7 km and 6 by 7.4 km (Figure 2b). The variable grid allows the use of 
more (less) points in regions of large (small) gradients. 

The 21 sigma levels used are shown in Figure 3 and Table 3. The sigma values 
range from zero at the surface to minus one at the bottom with the vertical grid spacing 
proportional to the ocean depth. The vertical resolution has been chosen to be higher near 
the surface and the bottom in order to resolve both the Surface Boundary Layer (SBL) 
and the Bottom Boundary Layer (BBL), which are important in coastal regions. To 
eliminate the time constraints for the vertical grid related to the higher resolution near 
surface, bottom and shallow waters, an implicit vertical time differencing scheme is used. 

The prognostic variables of the model are potential temperature, salinity, density, 
the three components of velocity, surface elevation, turbulence kinetic energy and length 
scale. The model has a split time step for the external and internal modes. The external 
mode solves the equations for the vertically integrated momentum equations. It also 
provides the sea surface and barotropic velocity components, and has a time step of 6 
seconds. The internal mode solves the complete 3D equations and has a time step of 180 
seconds. 
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A Smagorinsky formulation (Smagorinsky et al., 1965) is used for the horizontal 
diffusion in which the horizontal viscosity coefficients depend on the grid size, the 
velocity gradients and a coefficient. In this study a value of 0.2 was assigned to this 
coefficient, consistent with other POM studies (e.g, Ezer and Mellor, 1997). 

C. INITIALIZATION, FORCING AND BOUNDARY CONDITIONS 

The model was initialized with annual temperature and salinity values obtained 

from Levitus and Boyer (1994) and Levitus et al. (1994). Due to the short time of the 

model runs (less than 60 days), zero salinity and temperature fluxes have been prescribed 

at the ocean surface. The climatological surface temperature (Figure 4a) shows a decrease 

in temperature from the southwestern to the northeastern comer of the domain. The 

gradient increases to the north (which will be shown to be important for the formation of 

the Iberian Current). A cross-section for salinity at 36°N (Figure 4d). shows the 

Mediterranean Water signature at -1200 m depth with salinity values of -36.1 psu. In the 

upper 300 m the North Atlantic Ocean waters are found, with the characteristic high 

salinity values. A cross-section for temperature at the same location (Figure 4e) shows 

below -1000 m depth a downward sloping of the isotherms approaching the coast which 

is consistent with the presence of warm Mediterranean Waters and of a poleward flow. In 

a cross-section near ~41.5°N (northern boundary) the signature of the Mediterranean 

Water is found at -1000 m depth with a much less zonal extent but with the same value 

for the maximum salinity (Figure 4f). The waters at 41.5°N are generally -0.1 to 0.2 psu 

less haline than the waters at 36°N. The temperatures for the 41.5°N cross-section (Figure 
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4g) still show a strong Mediterranean Water influence but now at -1000 m depth, 
consistent with the high salinity values found at this location. A cross-section of salinity 
at the southern boundary (Figure 4b) shows the Mediterranean signature at -1200 m 
depth with salinity values of - 35.7 psu, much less than the ones found on the northern 
boundary. The temperature signature (Figure 4c) is also less strong than in the northern 
boundary case (the vertcal gradient of temperature near the coast is much higher in the 
southern boundary case). 

The model was forced from rest with the annual ECMWF wind fields, which were 
interpolated for the model grid. As expected, the wind stress is stronger in the southern 
region of the model domain and weaker off Iberia and in the Gulf of Cadiz (Figure 5). 

Correct specification of the open boundary conditions (BC) is very important to 
achieve realistic results, with no reflections, clamping, spurious currents or numerical 
alteration of the total volume of water in the model. The problem is that there is not a 
general criteria that can give the answer to what boundary conditions are the best for a 
specific model or study. For models with a free surface, such as used here, one of the 
important criteria is that the BCs should be transparent to the waves. In this model, a 
gradient boundary condition (Chapman, 1985), which allows geostrophic flow normal to 
the boundary, worked best for the elevation. For the baroclinic velocity components 
normal to the boundary, an explicit wave radiation scheme based on the Sommerfeld 
radiation condition was used. For inflow situations, the model was forced with annual 
temperature and salinity values obtained from Boyer (1994) and Levitus et al. (1994), 
while in outflow situations an advection scheme was used. 
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For the barotropic velocity components, a Flather radiation plus Roed local 
solution (FRO) was used. Palma and Matano (2000) showed good results with the FRO 
during BC tests to determine the BCs response to an alongshelf wind stress. Matano 
(1998) also showed that the FRO BC had good reflection properties and results in a test 
that determined the BC response to the combined action of wind forcing and wave 
radiation. His tests were executed with the barotropic version of POM and compared with 
benchmark results (no boundary conditions). 
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III. RESULTS FROM MODEL SIMULATIONS 

A. EXPERIMENT 1 - PRESSURE GRADIENT FORCE ERROR 

In Experiment 1 (see Table 1), the model was initialized with the horizontally 
averaged annual climatological temperatures and salinities. A realistic coastline and 
realistic topography were used, and there was no wind or thermohaline forcing. 

With the horizontal averages of the climatology and no forcing, we should expect 
that nothing will happen, i.e., the initial state of rest should be mantained with time. Due 
to pressure gradient force errors, however, this will not be the case and there will be 
resultant velocities as the result of these errors. 

Velocity errors induced by the pressure gradient force are unavoidable in 3D 
sigma coordinate models. There are two types of sigma coordinate errors, the sigma error 
of the first kind (SEFK) and of the second kind, as defined by Mellor (1998). The first 
one goes to zero prognostically by advecting the density field to a new state of 
equilibrium. The second one, a vorticity error, is the most important because it does not 
vanish with time, and is present in both 2D and 3D cases. 

There are several techniques to reduce the pressure gradient errors: 

1 - Smoothing the topography can reduce both SEFK and SESK. In particular, the 
slope parameter should not be greater than 0.2 (Mellor, 1998). Greater values of 
this parameter can induce currents over 1 m/s. 

2 - Using the highest possible resolution can reduce the errors, since, the pressure 
gradient error decreases with the square of the horizontal and vertical grid size 
(Mellor et al., 1994). 
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3 - Subtracting the horizontally averaged density before the computation of the 
baroclinic integral reduces the SESK (Mellor, 1998). 

4 - Using a curvilinear grid that follows the bathymetry reduces the SESK 
(Mellor, 1998). 

In this study the first three techniques were used. The last technique, the use of a 
curvilinear grid, could not be used, because the geography of the Gulf of Cadiz would 
have given rise to singularity points. 

The results of Experiment 1 are shown in Figures 6 and 7. In Figure 6 the velocity 
has been plotted versus time. The velocity error is shown to increase faster in the first 3 
days, and become almost constant after day 7. The size of the error reduces to less than 
0.5 cm/s by day 10. To show where the errors are present in the model domain, the 
velocity field at sigma level 1 is shown in Figure 7. As expected, maximum velocities of 
~0.5 cm/s are found within -30 km from the coast where the slope parameter is the 
largest. Experiment 1 has shown that with the use of the three techniques, the pressure 
gradient error has been considerably reduced. Before the use of these techniques, model 
runs showed pressure gradient errors of - 1 m/s in the coastal regions (not shown). 


B. EXPERIMENT 2 - WIND FORCING ON A FLAT BOTTOM 

In Experiment 2 (see Table 1), the model was initialized with the horizontally 
averaged annual climatological temperatures and salinities. A realistic coastline and flat 
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bottom (constant depth of 2500 meters) were used, and the model was forced with annual 
climatological winds. 

By day 10, a quasi-steady state for the kinetic energy has been achieved (not 
shown). The strong winds at the southward end of the model have caused cooler 
temperatures associated with strong upwelling in the coastal region south of Cape 
Beddouzza (e.g., Figure 8). Throughout the model domain, there are predominantly 
southward currents except in a zonal oceanic region between 32.5 N and 34.5 N where 
there is a slight northward component due to the Ekman transport (e.g, see wind stress. 
Figure 5). The stronger currents tend to be found near coastline features, e.g., off Cabo da 
Roca and off Cape Ghir, with maximum speeds of ~50 cm/s in the coastal region between 
Cape Beddouzza and Cape Ghir (e.g.. Figure 9a). Between 1000 m depth (Figure 9b) and 
the bottom (Figure 9c), the flow is predominantly barotropic throughout the domain (e.g., 
compare Figures 9b and 9c). The coastal current at these depths is poleward, in the 
opposite direction of the surface flow. 

By day 20, the offshore extent of the coastal upwelled waters has increased, 
particularly off the coastal area south of Cape Ghir where the offshore extent is -120 km. 
The coldest waters are found south of Cape Ghir and Cabo da Roca (e.g.. Figure 10). Off 
the Iberian coast and south of Cape Beddouzza, the waters are about 1°C colder than at 
day 10. The maximum surface velocities are found in the coastal region near Cape 
Bedouzza (not shown). 

By day 30, the temperature fields (Figure 11) show that the upwelled waters have 
propagated farther westward. The offshore extension of the 17°C isotherm is at -300 km 
off Cape Ghir and at -60 km off Iberia. The minimum offshore extension is in the Gulf of 
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Cadiz region. The surface temperatures near the coast are very similar to the ones at day 
10 . A filament is present off Cape Ghir with an offshore extension of ~55 km (e.g., 
Figure 7). Figure 12 shows the presence of disorganized flow near this area, which 
means that some of the mean kinetic energy has been converted to eddy kinetic energy. 

By day 40, there is evidence of filament activity off Cabo da Roca, the southwest 
tip of Iberia, and off Cape Ghir (Figure 13). The development of mesoscale features is 
also evident in almost all of the coastal domain with the more developed features near 
coastline irregularities, i.e., Cabo da Roca, Cabo de Sao Vicente, Cape Ghir and Cape 
Beddouzza. The development of some smaller features in the region of the Gulf of Cadiz 
is also discernible (Figure 14). 

Also apparent in the model results is the development of two anticyclonic eddies 
at -1200 m depth (e.g., Figure 15), one off Cabo de Sao Vicente and another one off 
Cabo da Roca. A cross-section of salinity at 37.42°N, near the southern anticyclonic eddy 
(e.g.. Figure 16) shows a salty core of 35.87 psu at 1100 m depth, which is associated 
with the signature of Mediterranean water (present in the initial annual climatology). It 
has been suggested (e.g., Kase et al., 1989) that the Meddies are generated by the basic 
instability of the equatorward coastal jet and the poleward undercurrent. The deep origin, 
salty signature, depth and anticyclonic rotation of the eddy west of Cabo Sao Vicente is 
consistent with observations of Meddies in this region (e.g., Richardson and Tychensky, 
1998). The Meddy observed off Cabo da Roca is consistent with the results of Kase et al., 
(1989), who observed eddies in the Mediterranean Outflow off Iberia. 

The continuing westward propagation of upwelled water with time, as seen in 
Figures 8, 10, 11 and 13, is consistent with the results of McCreary et al. (1987) and 
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Batteen et al. (1989). They showed that, due to the beta effect, the surface coastal jet does 
not necessarily have to be confined to within a Rossby radius of deformation of the coast. 

The results in this experiment are also consistent with the ones obtained by 
Batteen et al. (2000) with a primitive equation model with biharmonic diffusivity, based 
on the hydrostatic, Boussinesq and rigid lid aproximations. In particular they were able to 
reproduce a realistic equatorward current, a poleward undercurrent, Meddies and other 
mesoscale features (eddies, meanders and filaments). 

C. EXPERIMENT 3 - ROLE OF TOPOGRAPHY 

In Experiment 3 (see Table 1), the model was initialized with the horizontally 
averaged annual climatological temperatures and salinities. A realistic coastline and 
bottom topography were used along with forcing by the annual climatological winds. 

By day 10, the velocities in the coastal areas in the bottom topography 
Experiment 3 (e.g., Figure 17a) are roughly twice the ones in the flat bottom Experiment 
2 (e.g.. Figure 9a). The velocities in Experiment 3 also have the highest magnitudes at 
~15 km offshore (near the shelf break) instead of at the coast as in the flat bottom 
Experiment 2. As in the flat bottom experiment, the strongest currents are found in the 
coastal region south of Cape Beddouzza. The temperatures at the surface are also very 
similar for the two experiments (e.g., compare Figures 8 and 18), both in the location of 
the strongest upwelling regions (south of Cape Beddouzza) and in the horizontal extent 
of the upwelled waters. One of the main differences between the experiments occurs off 
Figueira da Foz where the upwelled waters are present farther from the coast in 
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Experiment 3. The main reason for this phenomena is the location of a topographic 
feature at this location. Below the surface, at intermediate sigma level 10, the coastal 
velocities are still equatorward but are about 50% of the ones at surface (Figure 17b) 
while near the bottom the equatorward velocities decrease to 10% of the surface values 
(Figure 17c). 

By day 20, a comparison of the 17°C isotherm in both experiments (e.g, see 
Figures 10 and 19) shows that the extent of the upwelled waters is much less than in the 
flat bottom case (Figure 10) particularly off Cape Ghir. This can be explained by the 
presence of the bottom topography which traps the flow and opposes the tendency for 
westward propagation due to the planetary beta effect. In Experiment 3 the surface 
temperatures are also colder than in Experiment 2 in the region of Gulf of Cadiz where 
the winds are not as strong. 

The plot of total kinetic energy for Experiment 3 (e.g. Figure 20a) shows a steady 
and rapid increase of the kinetic energy for the first 25 days. During this period the flow 
develops currents but no mesoscale features. Between days 25 and 35, the kinetic energy 
decreases to about 1/3 of the maximum value. This loss consistent with (R0ed, 1999) can 
be explained by a transfer of kinetic energy to available gravitational energy and then to 
eddy gravitational energy and eddy kinetic energy (e.g., Figure 20b). This explanation is 
supported by the development of eddies at ~ day 30. The eddies develop at Cabo da 
Roca, Cabo Sao Vicente and off Figueira da Foz (e.g, Figure 21). The surface 
temperature (e.g, Figure 22) is similar to experiment 2 with the waters being ~1°C colder 
in the coastal area south of Cape Ghir. Off Cape Ghir the extent of the upwelled waters is 

also smaller. Note that the filament present in this region in Experiment 2 is now absent 
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in Experiment 3. This shows that the trapping of the currents due to topography is much 
stronger than the combined effect of coastline irregularities and the planetary beta effect. 
Note the presence of relatively strong eddies near the capes off Iberia, which were much 
weaker in Experiment 2. This shows that the bottom topography plays an important role 
in their intensification. 

Unlike Experiment 2 which had meanders throughout the coastal domain by day 
40 (e.g see Figure 14), the only evidence of mesoscale phenomena in this experiment are 
two eddies, one off Cabo da Roca and the other off Figueira da Foz (e.g. Figures 23a and 
23b). The eddy off Cabo da Roca is centered at ~38.82°N, 10.44°W, with a diameter of 
70 km and a core temperature of 14.4°C. Off Cape Ghir and off the coast of Portugal 
between Cabo da Roca and Cabo Sao Vicente, the flow is mantained much closer to 
shore in the experiment with topography (e.g., compare Figures 13 and 23a). The 
opposite situation is shown in the Gulf of Cadiz where the slope of topography is much 

less than in the rest of the coastal domain. 

The poleward undercurrent is much weaker in Experiment 3 then in Experiment 
2. For example, at day 30, it is less than lcm/s (Figure 24a) at 1200 meters depth in 
Experiment 3, while in Experiment 2 the current was 6 to 7 cm/s at the same depth (see 
Figure 24b). The core of the undercurrent in Experiment 3 is also much deeper. For 
example, at day 30 the core is at -2300 m depth while in Experiment 2 it is at -1200 m 

depth. 

Meddy formation also occurs much later in Experiment 3 than in Experiment 2. 
Instead of having Meddies form by day 30 as in Experiment 2, Meddies do not form until 
- day 60. The Meddies form off Cabo da Roca and off Figueira da Foz (Figure 25); 
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however, unlike Experiment 2, no Meddy forms off Cabo Sao Vicente. As Experiment 2 
showed, the sharpness of the coastline is likely necessary to the formation of Meddies off 
Cabo Sao Vicente. In Experiment 2 the Meddies also appeared to be generated due to the 
curvature of the Mediterranean outflow back on itself, consistent with the theory 
suggested by Pichevin and Nof (1996). It is speculated that the smoothing of topography 
also smoothed the irregularities of the coastline at 1200 depth. 

D. EXPERIMENT 4 - THE ROLE OF BOTTOM TOPOGRAPHY WITH FULL 

CLIMATOLOGY 


In Experiment 4 (see Table 1), the model was initialized with the full annual 
climatological temperatures and salinities. A realistic coastline and bottom topography 
were used along with forcing by the annual climatological winds. 

By day 10, the currents are very similar to the ones found on Experiment 3 (i.e., 
compare Figures 26 and 17a), except ~80 km offshore of Iberia where there is poleward 
flow, and warmer temperatures off Cabo Sao Vicente, which are advected poleward (see 
the 18°C isotherm. Figure 27) . Since this is a region of weak wind stress, the Ekman 
transport is not able to overcome the geostrophic adjustment of density due to the 
thermohaline gradients from the full temperature and salinity fields. Even though the 
Gulf of Cadiz is also an area of weak wind stress, this phenomena is not present due to 
the fact that the horizontal variation of temperature (density) is much smaller (Figure 4a) 
compared to the Ekman transport. The difference in temperatures in coastal areas 
between Experiments 3 and 4 (Figures 18 and 27, respectively) reflect the different 
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climatologies used. To the south, i.e., off Cape Ghir, the temperatures are similar in both 
experiments while in the northern area, i.e., off Figueira da Foz, Experiment 4 
temperatures are ~1.5°C colder than those of Experiment 3. The temperature differences 
between the waters in upwelling regions and adjacent waters for both experiments are 
approximately constant. 

Again in day 20, similar patterns are shown for the velocities for the two 
experiments (Figures 28 and 29). The strongest currents are found in the coastal area 
south of ~34°N except for ~80 km offshore of Iberia, where the currents have a poleward 
component in Experiment 4. A comparison of cross-sections of the meridional velocity 
component for Experiments 3 and 4 (e.g., Figures 30 and 31), at 32.8°N shows a much 
more developed poleward undercurrent in Experiment 4 with a well defined core at -800 
m depth (Figure 30) and velocities 20% stronger than in experiment 3 (Figure 31). In 
Experiment 4 the equatorward current is narrower and reaches -1400 m depth, compared 
to -1000 m depth in Experiment 3. 

A comparison of cross-sections of meridional velocity off Cabo Roca for 

Experiments 3 and 4 shows the presence of the Iberian current in Experiment 4 (Figure 

32) but not in Experiment 3 (Figure 33). The Iberian Current is a seasonal shoaling of the 

poleward undercurrent off Iberia (Haynes and Barton, 1990), usually seen in winter time, 

when the seasonal wind stress is also weaker off Iberia. This shows that the horizontal 

variability of density is important for the formation of the Iberian Current. At this 

location, the undercurrent in Experiment 4 has its core at ~800m with a maximum value 

of 27 cm/s, while the equatorward current is trapped within the first 400 m depth and has 

an offshore extent of -70 km (Figure 32). In Experiment 3 the undercurrent is much 
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weaker with a magnitude of -6 cm/s, does not have a well defined core and is found 
much deeper. Above the undercurrent, the equatorward current is 25% stronger than in 
the previous case, reaches 1500 m depth and is -140 km wide. As a result, a much 
stronger along shelf variability is present in Experiment 4 than in Experiment 3. 
Consistent with Batteen et al., 2000, and Batteen and Murray (2000), the Iberian Current 
is seen as far south as Cabo Sao Vicente 37.4°N (e.g.. Figure 34). 

By day 30, as in Experiment 3, there is a 20% decrease in the coastal surface 
velocities compared to those at day 20 (compare Figures 35 and 28), which could be due 
to two factors: first, the geostrophic adjustment produces velocities in opposition (but 
smaller in magnitude) to the main flow, resulting in a reduction of the speeds; second, 
the horizontal variability of density contributes to the enhancement of the poleward 
alongshore variability, giving rise to a greater enhancement and surfacing of the 
undercurrent, which subsequently reduces the flow surface speeds. 

A cross section at 32.8°N (Figure 36) shows a strong undercurrent with a speed of 
-20 cm/s with a core depth at -500 m depth for Experiment 4, a core of similar 
magnitude but with a much deeper location is shown for Experiment 3, at 1000 m depth 
(Figure 37). Off Iberia, the horizontally opposing flows of the equatorward coastal and 
the poleward Iberian current produced horizontal shear that could decrease the overall 
speeds (Figure 35). A cross-section at 32.8°N for Experiment 3 at day 40 (Figure 38) 
shows that the undercurrent has a shallower core at -750 m depth and a weaker (from 
69cm/s to 55 cm/s, a reduction of 20%) surface equatorward current. 

The meanders present off Figueira da Foz and Cabo da Roca in Experiment 4 are 
consistent with those in Experiment 3 . The extent of the upwelled waters is very similar 
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in Experiments 3 and 4 (Figures 62 and 19) in the coastal regions south of Cabo Sao 
Vicente; however, to the north there are differences in extent as well in value due mainly 

to differences between the two climatologies. 

By day 40, surface velocities are similar to the ones at day 30, except for the 
stronger surface poleward flow off Iberia and north of Figueira Foz (Figures 40 and 35). 
By day 60, a Meddy has formed off Cabo da Roca; however, as in Experiment 3, no 
Meddy forms off Cabo Sao (Figure 41). As in Experiment 3, the Meddy forms at ~ day 
60, compared to ~ day 40 in Experiment 2. 

Surface temperature and velocity fields at day 60 are shown for Experiment 4 
(Figure 42), Experiment 3 (Figure 43) and Experiment 2 (Figure 44). While all three 
experiments show similar features, i.e., surface equatorward flow, upwelling, meanders 
and eddies, notable differences are apparent. A comparison of the flat bottom (Figure 44) 
and topography (Figure 43) with the same climatology shows the effect of planetary beta 
acting alone in the flat bottom Experiment 2 (Figure 44) versus competing with the 
topographic beta effect in Experiment 3. The effect of the full climatology is evident in 
Experiment 4 (Figure 42), which shows the thermal gradient from north to south, the 
tightening of the currents near the coast and slightly weaker surface currents due to the 
effect of thermohaline gradients. 
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IV. SUMMARY 


The objective of this study was to investigate the roles of wind forcing, bottom 
topography and thermohaline gradients on classical as well as unique features in the 
northern Canary Current System. Toward this end, four numerical experiments were run, 
all on a beta-plane, with a sigma coordinate numerical model, i.e., the Princeton Ocean 
Model. The first experiment investigated the pressure gradient force error. The second 
experiment investigated the effect of annual wind forcing on a flat bottom. The third 
experiment investigated the additional effect of topography. The fourth experiment was 
used to determine the additional role of the full annual climatology. 

Experiment 1, used to evaluate the pressure force gradient error, showed that 
velocity errors inherent in three dimensional sigma coordinate models could be 
sucessfully reduced from ~1 m/s to ~0.5cm/s using three techniques: smoothing the 
topography, using the highest possible resolution, and subtracting the area averaged 
density before the computation of the baroclinic integral. The results showed that the 
highest velocities (~0.5cm/s) were concentrated near the coast where the values of the 
slope parameter were the highest. 

Experiment 2 produced an offshore surface equatorward meandering jet, realistic 
surface and subsurface poleward currents, upwelling, meanders, eddies and filaments. In 
addition, these experiments depicted unique NCCS features, including the geographical 
separation of the Gulf of Cadiz region from the west coast upwelling regimes, poleward 
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spreading of the Mediterranean Outflow, and the development and propagation of 
Meddies from the Cabo Sao Vicente and Cabo da Roca regions. 

A comparison between Experiments 2 and 3 showed that the bottom topography 
plays an important role in trapping and intensifying the equatorward current near the 
coast. The poleward undercurrent also weaker and had a deeper core depth in Experiment 
3. Stronger eddies occurred off Cabo da Roca and Cape Ghir in Experiment 3. No 
formation of Meddies off Cabo Sao Vicente in Experiment 3 occurred. It was speculated 
that since Meddies did form off Cabo Sao Vicente in Experiment 2, consistent with 
observations, that the sharpness of the coastline is likely to be necessary to the 
formation of Meddies. 

In Experiment 4, the additional effect of the full annual climatology produced an 
offshore surface equatorward meandering jet, realistic surface and subsurface currents, 
upwelling, meanders and filaments. The effect of the full climatology produces the 
tightening of the currents near the coast and slightly weaker currents due to the effect of 
thermohaline gradients. As in Experiment 3, there was no development of eddies of Cabo 
Sao Vicente but only off Cabo da Roca. Only Experiment 4 showed the Iberian Current 
off the Portugal coast showing that thermohaline gradients are essential to the formation 
of this current. 
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Figure 1. The model domain for the northern Canary Current System (NCCS) is bounded 
by 31.5°N to 41.5°N, 6°W to 16.5°W. Geographical locations (with abbreviated names 
used in the following figures) and prominent features are labeled. The model domain has 
a closed boundary along the entire coast and three open boundaries. 
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Figure 2a. Original topography (from Sandwell D.T. and W.F. Smith, 1996), with a 
resolution of 2 minutes (i.e., 1/30 of a degree). 
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Figure 2b. Grid lines with every fifth grid line plotted. 


29 






Figure 2c. Smoothed topography obtained after applying a linear two-dimensional low- 
pass filter and reassigning 2500 m depth to depths greater than 2500 m. 


30 











-16 -15 -14 -13 -12 -11 -10 -9 -8 -7 -6 

Longitude 


Figure 4b. Cross-section at 31.5°N of Levitus annual climatological salinity. 
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Figure 4c. Cross-section at 31.5°N of Levitus annual climatological temperature. 
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Figure 4d. Cross-section at 36°N of Levitus annual climatological salinity. 








Lat: 36N Temperature 



Figure 4e. Cross-section at 36°N of Levitus annual climatological temperature. 
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Figure 4g. Cross-section at 41.5°N of Levitus annual climatological temperature. 
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Figure 5. Wind stress in Pascal calculated from annual climatological ECMWF winds 
obtained from Trenberth et al. (1990). 
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Figure 9a. Surface velocity vectors for Experiment 2 on day 10. 
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Figure 9b. Velocity vectors at sigma level 10 for Experiment 2 on day 10. 
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Exp: 2 Day: 10 Bottom Velocity 



Figure 9c. Velocity vectors at sigma level 20 for Experiment 2 on day 10. 
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Experiment 2 on day 20. Contour interval is 
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Figure 12. Surface velocity vectors for Experiment 2 on day 30. 
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Exp: 2 Day: 40 Surface Velocity 



Figure 14. Surface velocity vectors for Experiment 2 on day 40. 
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Figure 15. Velocity vectors at 1200 m depth for Experiment 2 on day 40. 
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Figure 17a. Surface velocity vectors for Experiment 3 on day 10 
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Figure 17b. Velocity vectors at sigma level 10 for Experiment 3 on day 10. 
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Exp: 3 Day: io Surface Temperature 



Figure 18. Surface temperature contours for Experiment 3 on day 10. Contour interval is 
1 °C. 
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Figure 20b. Energy conversions in percentage of the wind work from Roed (2000). 
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Exp: 3 Day: 30 Surface Velocity 



Figure 21. Surface velocity vectors for Experiment 3 on day 30. 
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Figure 23a. Surface temperature contours for Experiment 3 on day 40. Contour interval 
is 1°C. 
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Figure 23b. Surface velocity vectors for Experiment 3 on day 40. 
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Figure 24a. Cross-section of meridional velocity (v) at 37.4°N for Experiment 3 on day 
30. Equatorward (poleward) flow is denoted by dashed (solid) lines with contour 
intervals of 5 cm/s (2cm/s). 
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Figure 27. Surface temperature contours for Experiment 4 on day 10. Contour interval is 
1°C. 
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Exp: 3 Day: 20 Surface Velocity 



Figure 29. Surface velocity vectors for Experiment 3 on day 20. 
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Figure 30. Cross-section of meridional velocity (v) at 32.8°N for Experiment 4 on day 
20. Equatorward (poleward) flow is denoted by dashed (solid) lines with contour 
intervals of 5 cm/s (2cm/s). 
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Figure 31. Cross-section of meridional velocity (v) at 32.8°N for Experiment 3 on day 
20. Equatorward (poleward) flow is denoted by dashed (solid) lines with contour 
intervals of 5 cm/s (2cm/s). 
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Figure 32. Cross-section of meridional velocity (v) at 39°N for Experiment 4 on day 20. 
Equatorward (poleward) flow is denoted by dashed (solid) lines with contour intervals of 
5 cm/s (2cm/s). 
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Figure 33. Cross-section of meridional velocity (v) at 39°N for Experiment 3 on day 20. 
Equatorward (poleward) flow is denoted by dashed (solid) lines with contour intervals of 
5 cm/s (2cm/s). 
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Exp: 4 Day: 30 Surface Velocity 



Figure 35. Surface velocity vectors for Experiment 4 on day 30. 
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Exp: 4 Day: 30 Lat: 32.8N Velocity 
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Cross-section of meridional velocity (v) at 32.8°N for Experiment 4 on day 30. 
Equatorward (poleward) flow is denoted by dashed (solid) lines with contour intervals of 
5 cm/s (2cm/s). 













Longitude 


Figure 37. Cross-section of meridional velocity (v) at 32.8°N for Experiment 3 on day 
30. Equatorward (poleward) flow is denoted by dashed (solid) lines with contour 
intervals of 5 cm/s (2cm/s). 
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Figure 38. Cross-section of meridional velocity (v) at 32.8°N for Experiment 3 on day 
40. Equatorward (poleward) flow is denoted by dashed (solid) lines with contour 
intervals of 5 cm/s (2cm/s). 

















Exp: 4 Day: 30 Surface Temperature 



Figure 39. Surface temperature contours for Experiment 4 on day 30. Contour interval is 
1°C. 


80 






Exp: 4 Day: 40 


Surface Velocity 



Longitude 


Figure 40. Surface velocity vectors for Experiment 4 on day 40. 
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Figure 41. Surface velocity vectors (arrows) and temperatures (in color) for Experiment 
4 on day 60. 
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Figure 43. Surface velocity vectors (arrows) and temperatures (in color) for Experiment 
3 on day 60. 
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Figure 43. Surface velocity vectors (arrows) and temperatures (in color) for Experiment 
2 on day 60. 
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Depth (m) 

Level 

Depth (m) 

Level 

Depth (m) 

T 

0 

12 

300 

23 

1400 

2 

10 

13 

400 

24 

1500 

3 

20 

14 

500 

25 

1750 

4 

30 

15 

600 

26 


5 

50 

16 

700 

27 

2500 

6 

75 

17 

800 

28 

HESSSm 

7 


18 

900 

29 

■ 

8 

125 

19 

1000 

30 

4000 

9 

150 

20 

1100 

31 

4500 

10 

200 

21 

1200 

32 

5000 

11 

250 

22 

1300 

33 

5500 


Table 2. Vertical levels and depths used by Levitus and Boyer (1994) 

and Levitus et al. (1994) 

I 

I 


Level 

Sigma Value 

Level 

Sigma Value 

1 

0 

12 

-0.61538 

2 

-0.00961 

13 

-0.69231 

3 

-0.01923 

14 

-0.76923 

4 

-0.03846 

15 

-0.84615 

5 

-0.07692 

16 

-0.92308 

6 

-0.15385 

17 

-0.96154 

7 

-0.23077 

18 

-0.98077 

8 

-0.30769 

19 

-0.99038 

9 

-0.38462 

20 

-0.99519 

10 

-0.46154 

21 

-1.0 

11 

-0.53846 




Table 3. Values of sigma levels 
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